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^ ' Abstract 

£N| . Cylindrical algebraic decomposition (CAD) is an important tool for the investigation of semi-algebraic 
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sets. Originally introduced by Collins in the 1970s for use in quantifier elimination it has since found 
^ ] numerous applications within algebraic geometry and beyond. Following from his original work in 1988, 

McCallum presented an improved algorithm, CADW, which offered a huge increase in the practical utility 
of CAD. 

In 2009 a team based at the University of Western Ontario presented a new and quite separate 
algorithm for CAD, which was implemented and included in the computer algebra system Maple. As 
part of a wider project at Bath investigating CAD and its applications, Collins and McCallum's CAD 
, algorithms have been implemented in Maple. This report details these implementations and compares 

X^/y ' them to Qepcad and the Ontario algorithm. 

The implementations were originally undertaken to facilitate research into the connections between 
the algorithms. However, the ability of the code to guarantee order-invariant output has led to its use 
in new research on CADs which are minimal for certain problems. In addition, the implementation 
described here is of interest as the only full implementation of CADW, (since Qepcad does not currently 
make use of McCallum's delineating polynomials), and hence can solve problems not admissible to other 
CAD implementations. 

O ; This work is supported by EPSRC grant EP/J003247/1. 

csj ' 1 Introduction 

O ' 

£f~) • Cylindrical algebraic decomposition (CAD) was first announced in 1973 by Collins. A CAD is a decomposition 

of R™ into cells, constructed with respect to a set of input polynomials in n variables. Each cell can be 
described as a semi-algebraic set and the cells are cylindrically arranged, meaning the projection of any two 
cells is either equal or disjoint. Usually, the CAD produced is such that each polynomial is sign-invariant 
within each cell, allowing for the solution of many problems defined by the polynomials. Collins provided 
the first algorithm to compute a CAD jTTJ [T] , developed as a tool for quantifier elimination in real closed 
fields. Since their discovery they have found numerous applications ranging from robot-motion planning |20j 
to simplification technology [H [14j etc.] 

Collins' algorithm has two phases. The first, projection, applies a projection operator repeatedly to a 
set of polynomials, each time producing another set of polynomials in one fewer variables. Together these 
sets contain the projection polynomials. The second phase, lifting, then builds the CAD incrementally from 
these polynomials. First R is decomposed into cells which are points and intervals corresponding to the 
real roots of the univariate polynomials. Then R 2 is decomposed by repeating the process over each cell 
using the bivariate polynomials at a sample point of the cell. The output for each cell consists of sections 
of polynomials (where a polynomial vanishes) and sectors (the regions between these). Together these form 
the stack over the cell, and taking the union of these stacks gives the CAD of R 2 . This process is repeated 
until a CAD of R™ is produced. 

To conclude that the CAD of R" produced using sample points in this way is sign- invariant we need 
the key concept of delineability. A polynomial is delineable in a cell if the portion of its zero set in the cell 
consists of disjoint sections. A set of polynomials are delineable in a cell if each is delineable and further that 
the sections of different polynomials in the cell are either identical or disjoint. A projection operator is valid 
for use in the algorithms if over each cell of a sign- invariant CAD for projection polynomials in r variables, 
the polynomials in r + 1 variables are delineable. 

The output of a CAD algorithm depends on the ordering of the variables. In this paper we usually work 
with polynomials in Z[xi, . . . ,x n ] with the variables, x, listed in ascending order; (so we first project with 
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respect to x n and so on until we have univariate polynomials in x\). The main variable of a polynomial, 
mvar(/), is the greatest variable present with respect to the ordering. 

Since Collins published the original algorithm there has been much research into improvements. These 
include, but are not restricted to; the use of equational constraints [TH1 etc.], the notion of partial CADs [T31 
etc.] and [HJ etc.], ideas on cell adjacency and clustering [H etc.], ideas on preprocessing the input [13J etc] 
and heuristics to pick the best variable ordering |15[ etc.]. A summary of developments over the first twenty 
years of CAD was given by [12] . 

A key area of research has been in improving the projection operator used in the first phase of the 
algorithm. Collins original operator can produce CADs with many more cells than required for sign invariance 
of the polynomials. McCallum proved in |17j that for most problems a far simpler projection operator could 
be used. Then in [15] McCallum presented an algorithm, named CADW, detailing how and when the improved 
operator should be used. A particular feature of the algorithm is that the CAD produced is not just sign 
invariant with respect to the input polynomials, but order invariant. Further development of McCallum's 
algorithm and operator have been presented by Brown in [SJ [3] • 

McCallum's projection operator is implemented in Qepcad [7]; an interactive command-line program 
written in C using the SACLIB library of computer algebra functions. The name stands for quantifier elimi- 
nation by partial cylindrical algebraic decomposition but the software can also produce full cads and offers a 
wide range of options such as the use of equational constraints and 2d plots. 

In 2009 a new approach to producing CADs was presented in [TU]. Instead of projection and lifting this 
approach first computes a decomposition of complex space and then refines to real space, making use of 
the theory of triangular sets and regular chains [3J etc.]. The algorithm was implemented in Maple and is 
included in the RegularChains library [T^l etc.], distributed with Maple. 

This report gives details on a new implementation of CAD via projection and lifting, within Maple. 
The implementation was originally undertaken to facilitate research into the connections between the two 
approaches for computing CAD. However, its ability to guarantee order-invariant output has led the imple- 
mentation to utility as part of new algorithm in [5] offering decompositions which while not sign invariant 
for their input polynomials, arc minimal for certain problems based on them. In addition, the implementa- 
tion is of interest as it is the only full implementation of CADW, since Qepcad can not produce McCallum's 
delineating polynomials leading to unnecessary warnings about potential failure for certain examples. (We 
note that SyNRAC [M] also implements CAD via projection and lifting in Maple, but this implementation 
uses Collins' algorithm and so can not provide order- invariance for use in [5].) 

In Section [5] we summarize the theory of the Collins and McCallum projection operators algorithm and 
in Section [3J we describe our implementation in Maple. In Section 2] we describe how this implementation 
compares with both Qepcad and the algorithm from the theory of regular chains. This report should be 
accompanied with a file ProjectionCAD.mm containing the code implementing the algorithms in Maple. 
The file is designed to be read into Maple to define a library of commands, but can also be viewed as a text 
file. Details on installing and working with the code are given in Subsection 13. II 

2 Projection operators 

We summarize the projection operators of Collins and McCallum, (for more details see PQ and [T7]V 

Let / and g be real univariate polynomials. Denote by deg(/) the degree of /, by /' the first derivative 

of /, by red fc (/) the fcth reductum of /, by discr(/) the discriminant of /, by res(/, g) the resultant of / and 

g and by psc^- (f,g) the jth principal subrcsultant coefficients of / and g. 

Let F be a set of univariate polynomials and suppose /, g E F. Then denote by coeff(F), the set of all 

non-zero coefficients of all elements of F, by red(F) the set of all red fc (/) such that < k < deg(/), by 

psc(-F) the set of all psCj(f,g) such that / ^ g and < j < min(def (/), deg(g)), by psd(F) the set of all 

psCj(/, /') such that < j < deg(/) — 1, by discr(F) the set of all discr(/) and by res(F) the set of res(f,g). 

We can now define the operators 



Note that P(F) is a subset of PROJ(F). When dealing with multivariate polynomials we may think of them 
as univariate polynomials in the main variable and use the definitions as above. 

When Collins first devised CAD he showed that P may be used as the projection operator for problems in 
2 variables. However, P cannot always be used for problems in an arbitrary number of variables and so [11] 



PROJ(F) 
P(F) 



coeff(F) U psd(red(F)) U psc(red(F)), 
coeff(F) U discr(.F) U res(F). 
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specified the algorithm using PRO J instead. In [T7] McCallum proved that P may also be used for problems 
in 3 variables and then in |18| McCallum extended his results to explain when they hold for problems in 
more variables. 

Key to McCallum's proofs is the concept of order invariance. Recall that a polynomial has order of 
vanishing k at a point, if k is the smallest non- negative integer such that some partial derivative of the 
polynomial of order k docs not vanish at the point. A CAD is order invariant with respect to a set of 
polynomials if each polynomial has constant order of vanishing within each cell. 

Let F be a square free basis of polynomials in r variables. McCallum's key theorem showed that if P(F) 
is order invariant over a cell of R r_1 then each element of F either vanishes identically or is order invariant 
in every cell of the stack above. Hence the theorem can be applied recursively to prove P is valid, so long as 
the nullification of a polynomial over a lower-dimensional cell does not occur. 

To classify when the operator P may be used, McCallum defined the following condition. Let F be a set 
of polynomials in r variables and let prim(.F') and cont(-F) be respectively the set of primitive parts and the 
set of contents of F. Then F is well-oriented if cither r — 1 or both the following hold: 

(1) Every / G pvvaa(F) has a finite number of nullification points, i.e. if the main variable is Xk then 
f(a,Xk) = for at most a finite number of a G R fc_1 . 

(2) The set P(F) U cont(F) is well oriented. 

If condition (1) is changed to the stronger condition of no nullification points then the polynomials are said 
to be very well oriented. 

For well-oriented sets of polynomials the projection operator P can be used to generate CADs. This is 
shown using McCallum's theorem discussed above and, in the case where are a finite number of nullifying 
points, by replacing the nullified polynomial / with a delineating polynomial; a partial derivative of / 
whose roots include the Xk coordinate of every point where the order of / differs from the minimal order 
induced by a. In [5] Brown defined the minimal delineating polynomial as the case where these are the 
only roots. 

In [18] McCallum gave an algorithm CADW which aims to build a CAD using P. Before lifting over each 
cell the algorithm checks for nullification over the cell. If the dimension of the cell is zero then a delineating 
polynomial is used instead and otherwise the algorithm fails, declaring the polynomials not well-oriented. If 
the algorithm finishes then the CAD produced is order-invariant. 

Qepcad offers both McCallum's operator and Hong's modification of Collins' [13]. However, Qepcad 
does not implement the calculation of delineating polynomials and so can fail unnecessarily when using 
McCallum's operator in cases where the polynomials are well oriented but not very well-oriented. 

However, this situation is minimised since before declaring failure Qepcad ensures this is really necessary 
by running through some checks detailed in [8]. For example, it may be the case that nullification does not 
lead to order-invariance, (meaning the minimal delineating polynomial is a constant). Consider / = zy — x 
which is nullified by the point (x, y) = (0, 0) and so is well oriented but not very well oriented. However, for 
all values of z this polynomial has order 1 and so no delineating polynomial need be added, (and Qepcad will 
not fail). Other checks performed by Qepcad include; seeing if a delineating polynomial is already included 
in the projection set, checking whether the nullified polynomial is one that must actually be order-invariant 
to ascertain correctness, and checking whether it is the final lift in which case only sign-invariance is required. 

These checks rule out many cases where we may expect Qepcad to fail but the example constructed 
below demonstrates that there are still cases where delineating polynomials are really required to guarantee 
correct output using McCallum's operator. 

Example 2.1. To construct the example we need a polynomial in k variables nullified by a cell in k — 1 
variables. We must have k > 2 or the nullification will be avoided when considering the content and primitive 
part individually. A simple choice is f = zy — x 2 which is nullified over a cell with x — y = 0, which would 
certainly be produced from its coefficients belonging to the projection set. Further, this f has non- constant 
minimal delineating polynomial, z. 

However, if we were to just construct a CAD for f then the nullification would occur at the final lift and 
so we need not bother with a delineating polynomial. We must instead consider a polynomial p which produces 
f as a projection factor. If we let f be a coefficient then we need only make f order-invariant J6}/ so instead 
we choose a polynomial for which f is the discriminant: p = f + w 2 . 

Hence, when considering a CAD for p with respect to the variable ordering w > z > y > x, a delineating 
polynomial will certainly be required. This example was considered in the experiments detailed in Section ^ 
and using McCallum's operator a CAD with 73 cells was produced. 

1 Actually Qepcad does not necessarily fail, but it will produce a warning message indicating that the output may not be 
correct. 
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3 Implementation in Maple 

Algorithms to construct CADs with both Collins' and McCallum's projection operators have been imple- 
mented in Maple, with pseudo code for the main algorithms presented below and the actual code freely 
available from the author's website. The package is titled ProjectionCAD and is designed to complement 
the alternative approach of constructing CADs via regular chains which is already distributed with Maple. 
Indeed, the code described here makes use of tools from the RegularChains package and gives output in the 
same format. 

Algorithm Q] computes a set of projection polynomials (including the input polynomials) with respect to 
the chosen operator. Algorithm [2] will use those polynomials and the projection choice to build a CAD of W l 
over which they are sign invariant. Algorithm [3] is a composition of the two which can generate a variety of 
CADs via projection and lifting. 



Algorithm 1: CADProjection 



Input : • The input set of polynomials F C ]R[x„, . . . , xi], 

• A choice proj of cither Collins or McCallum. 
Output: • A set of projection polynomials P C R[x„, . . . , x{\. 

1 Set Pi to be the finest square-free basis for the primitive parts of F 

2 Set cord to be the set of contents of F. 

3 for j = 2,...,ndo 
Set Pi to be proj (Pi) U cord excluding any constant polynomials. 
Set Pi to be the finest square- free basis for the primitive parts of Pj. 
Reset cont to be the set of contents of Pi. 



r P Pi 
8 return P 



Both Collins and McCallum note in their work that for some examples simplifications could be made to the 
coefficients / rcducta included in their operators. For example, if a polynomial has a constant coefficient then 
this clearly does not need to be added to the projection polynomials. Further, no subsequent polynomials 
would need to be added since there will be no situations where these become leading coefficient. Some such 
elementary simplifications have been incorporated in the Maple implementation of Algorithm [TJ 

We note that Algorithm[2]diffcrs depending on the choice of projection operator in Algorithm[TJ due to the 
requirements to check that the input polynomials are well-oriented when using McCallum's operator. There 
is a further optional input in Algorithm [2] which can be used to request that the CAD outputted is not only 
sign- invariant but also order invariant. If not set to true then the algorithm does not check for nullification 
on the final lift since without this the CAD of R n is still guaranteed to be sign-invariant, sufficient for most 
applications, (but insufficient for the application in [5]). 

Algorithm |4] is a sub-algorithm which describes how the stack generation is implemented. To generate a 
stack the real roots of those projection polynomials with main variable Xi need to be calculated when the 
other variables are set according to a cell of dimension i — 1. New cells of dimension i are then defined with the 
variable Xi set in turn to be those roots and the intervals between. The roots may not be rational but algebraic 
numbers, and so care needs to be taken in implementing this. We make use of the RegularChains library 
and in particular and internal command for stack generation described in Section 5.2 of [10]. Algorithm |4] 
ensures that the polynomials passed to the RegularChains algorithm satisfy the assumptions that algorithm 
makes, namely that the polynomials are co-prime and square-free when evaluated on the cell, (separate above 
the cell in the language of regular chains). 

3.1 Working with the ProjectionCAD library in Maple 

This report should be accompanied by a file ProjectionCAD. mm containing the code implementing the algo- 
rithms in Maple. The file is designed to be read into Maple to define a library of commands, but can also 
be viewed as a text file. To read the package into Maple and make the commands available use 

> readC'ProjectionCAD .mm") : 

> with(ProjectionCAD) : 

Commands implementing Algorithms HH2] are now available, with the same names. For each there are two 
required arguments; a list of polynomials and a list of variables in descending order. There are also optional 
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Algorithm 2: CADLifting 



Input : • A set of polynomials P C K[x n , ■ • ■ , Xi] from CADProjection. 

• A choice proj of either Collins or McCallum. 

• A boolean fcad if the final CAD is to be order invariant. 
Output: A sign-invariant CAD of E™, also order-invariant if fcad = true. 

for i = 1 , . . . , n do 

Pi := {p G P such that the main variable of p is a;,}. 

Set Ci to be a CAD of R formed by the decomposition of the real line according to the real roots of 
the polynomials in P„ . 
for j = 2,...,ndo 

for each cell c G Q_i do 

if proj = McCallum then 

if i < n or fcad — true then 
set Qi to be the empty set. 
for each polynomial p G P n +i-i do 
if p = throughout c then 
if dim(c)=0 then 

add the minimal delineating polynomial to Qi if it is non-constant, 
else 

give warning message about potential failure. 

else 

add p to Qi 



4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 

15 
16 



17 

18 



19 
20 



else 

L Set Q t := Pi 

else 

L Set Q, = Pi 

Set S c := CADGenerateStack(c, Qi 
Set d := U c S c 
23 return C n 



Algorithm 3: CADFull 



Input : • The input set of polynomials F C R[xi, . . • , i n ]- 

• A choice proj of either Collins or McCallum. 

• A boolean fcad if the final CAD needs to be order invariant. 
Output: C, a CAD of E™. 



1 P := CADProj ection( F, proj) 

2 C := CADLift(P, proj, fcad) 

3 return C 



Algorithm 4: CADGencrateStack 

Input : • A cell c from a CAD of E i_1 , or if i = 0. 

• A set of projection polynomials P C R[x ra _i, . . . , x n ]. 
Output: S, a stack over c with respect to P. 

1 Set P to be a set of polynomials with the same zeros as P but which are coprime and square-free 
within the cell c. 

// By encoding the variables which are set to a point as a regular chain, and using commands from 
the RegularChains package [T5] . 

2 5:= RegularChains-GenerateStack(c, P) 

3 return S 
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keyword arguments; method to choose between Collins or McCallum, output to specify how the CAD is 
represented and f inalOI to specify that the outputted CAD be order invariant. Details on the progress of 
the algorithm can be accessed by setting the inf olevel for each command. 

Every cell in a cad of R™ is equipped with at least an index and a sample point. The index is an n-tuple 
of integers indicating the value or range of each variable using the real roots in increasing order computed for 
each stack. (Hence a cell whose index consists of only even integers defines a single point of R n while an index 
of only odd integers is a cell of full dimension n.) The sample points produced by Proj ectionCAD arc encoded 
as algebraic numbers given by regular chains and bounds isolating a single root. This is a consequence of 
using the internal RegularChains algorithm and is desired as it makes comparing CADs produced by the 
two approaches simple. 

The implementation in Maple can also produce cell representations (bounds or values for the variables 
using algebraic numbers) and if requested (with the output argument) can display the output intuitively 
using the piecewise construct in Maple. The number of cells in the CAD may be measured using the 
nops command on the output, (or the CADNumCellsInPiecewise command if the piecewise output format 
is used). 

A simple example of using the code is given below. The output is as displayed in a Maple worksheet, except 
that the sample points have been replaced by SP for brevity. 

> f := x~2+y~2-l: 

> cad := CADFull([f], vars, method=McCallum, output=piecewise) ; 

x < -1 



SP 



( SP 


V < 


I SP 


V = 


{SP 


< 


'sp 




SP 




< SP 


And ( 


SP 




SP 




(SP 


y < 


I SP 


y = 


{SP 


< 



-1 



y < —\/—x 2 + 1 
y = —\/—x 2 + 1 



y = +V—X 2 + 1 
v 7 — x 2 + 1 < y 



SP 



x = 1 

1 < X 



> CADNumCellsInPiecewise (cad) ; 



13 



4 Experimental results 

We have run experiments comparing three of the implementations of CAD. 

All the tests were run on a Linux desktop with Intel core i5 CPU (1.6GHz) and 8.0Gb total memory. The 
timings in Maple were recorded using the inbuilt time function while the timings for Qepcad are the total 
system time recorded at the end of each session. For the tests in Maple we use the timelimit function to 
limit each example to 1000 seconds of computation time. For the tests in Qepcad we run with the option 
+N500000000 + L200000, where the first option specifies the memory to be pre-allocated and the second the 
number of prime numbers to be used. 

The examples for the experiments were all taken from the CAD example repository described in |22j and 
stored at http://opus.bath.ac.uk/29503 and . Some of the problems here are quantified, however these 
experiments are designed to compare the CAD implementations only, and so were run on the unquantificd 
versions, (i.e. full sign- invariant CADs were computed for the polynomials involved in the problems). 
The three implementations compared were as follows. 



2 The Maple timelimit of 1000s does not apply to kernel operations so it is possible to have timings that are slightly higher, 
as is the case for the simplified Putnam example in Table [T] 
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PCAD We let PCAD denote the implementation of CAD via projection and lifting which is the topic of this 
paper. Run in Maple 16 through the command CADFull from the ProjectionCAD package developed 
at Bath. We run with the default settings; using McCallum projection and worrying about nullification 
at all except the final lift where this is ignored. 

TCAD We let TCAD denote the implementation of CAD via triangular sets. Run in Maple 16 with the 
in-built command from the RegularChains package, CylindricalAlgebraicDecompose. 

QCAD We let QCAD denote the implementation of CAD via projection and lifting given by Qepcad- 
B version 1.69. The code is given the full-cad option and default settings otherwise; McCallum 
projection and lifting with various improvements 6 ( 8 but no delineating polynomials. 

We note that there are other implementations of CAD available, including Mathematica, Redlog and SyN- 
RAC. 

The results of the experiments are displayed in Table [TJ An F indicates failure for theoretical reasons 
while a T/O indicates failure due to timeout in Maple. We note that with the setting described above, 
Qepcad never failed due to timeout or memory issues. 

Unless it fails for theoretical reasons, QCAD is usually the quickest. The examples where it was slower 
were generally small examples in which the time for QCAD was dominated by its initialisation. It is likely 
that Qepcad still represents the state of the art for computing full CADs at the moment. However, we 
make the cautionary note that the algorithm for TCAD described in [TU] is to be replaced by a significantly 
more efficient version, designed to reduce repeated computations and perform some steps in an incremental 
manner. This improvement is described in the preprint [5], but the implementation is still under development 
and not yet available. 

Next we note that TCAD never fails for theoretical reasons and is usually (but not always) quicker than 
PCAD. Comparing PCAD with QCAD we see that the cell counts are usually identical, as expected since 
they implement the same theoretical algorithm. There are occasions when QCAD will fail when PCAD 
does not because it does not implement McCallum's delineating polynomials, as in the example from Section 
[2] and the Quartic example. However, there are also examples where PCAD failed when QCAD did not, 
(the Whitney umbrella and the x-axis ellipse problem). These theoretical failures were avoided in QCAD 
since they correspond to cases where nullification does not contradict the theory underpinning the algorithm, 
which Qepcad checks for (as discussed in Section [5]). It is possible to instruct ProjectionCAD to continue 
the computation following a failure, which for these examples would lead to the same number of cells as 
QCAD. The other examples where the cell count differs are due to QCAD using a partial implementation of 
the simplified operator described in [6]. 

5 Summary 

We have described an implementation of CAD via projection and lifting in Maple. The implementation 
offers both McCallum's and Collins' operators with a variety of customizations available. The implementation 
was undertaken to facilitate research into the connections between the different approaches to CAD. The 
experimental results demonstrate that for most problems Qepcad is superior. Nevertheless, the approach 
described here offers significant utility for two main reasons. 

(1) It implements McCallum's delineating polynomials (actually Brown's minimal delineating polynomial) 
and thus can produce the only full implementation of McCallum's CADW. We have demonstrated that 
there are examples in which Qepcad will fail for which this implementation can succeed. 

(2) It can offer the user the choice of an order-invariant CAD in the final output. This is essential for the 
application in |5] which defines a new type of CAD, minimal for certain problems. 

Examples have been given where the output of CAD via Regular Chains in Maple is not order-invariant and 
it is not currently known how the algorithm could be modified to produce ordcr-invariance, (further research 
is currently being undertaken). 

It is likely that a small tweak order the hood in Qepcad could be made to offer order-invariant output 
at in (2). However, this would then greatly amplify the utility described in (1) since the majority of cases of 
nullification occur at the final lift where they can be safely ignored if only a sign- invariant CAD is required. 
For example, consider again / = zy — x 2 from Section [5] We noted there that / has non-constant minimal 
delineating polynomial and so to give an order-invariant CAD of / alone, this would need to be included. 



this includes the production of warnings that the output may not be correct 
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Running PCAD gives a CAD of 21 cells (as does TCAD and QCAD). One of these cells has x = 0, y = 0, x 
free, over which / is not order-invariant. However, if we run CADFull using the optional argument to specify 
an order-invariant output then this cell is split and the CAD outputtcd has 23 cells. 

If this implementation remains essential for applications such as [5] then some future improvements to 
the package arc likely. These should include the implementation of some of the improvements described in 
[5] to avoid failure unless it is absolutely necessary. Further, an implementation of the Brown-McCallum 
operator [6] is very desirable to offer the lowest possible cell counts. 
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M. England 



An implementation of CAD in Maple utilising McCallum projection 



Problem 


PCAD 


TCAD 


QCAD 


cells 


time 


cells 


time 


cells 


time 


Parametric parabola 


115 


0.5 


27 


0.1 


115 


1.6 


Whitney umbrella 


F 


- 


895 


3.4 


895 


3.8 


Cjuartic 


333 


3.9 


233 


2.7 


F 


- 


Sphere & catastrophe 


509 


6.2 


421 


4.140 


509 


3.7 


A ^ „ OA 

Arnon-84 


55 


0.3 


55 


0.157 


55 


3.7 


Arnon-84- 2 


41 


0.4 


41 


0.332 


41 


3.7 


Implicitization 


b 




one 

895 


5.228 


1 




Ball & cylinder 


365 


4.2 


365 


3.914 


365 


3.8 


Term rewrite system 


1099 


8.2 


1099 


7.768 


1099 


3.6 


Collins and Johnson 


3673 


50.2 


3673 


65.438 


3673 


3.8 


Lower bounds range 


F 




333 


1.550 


F 




X-axis ellipse problem 


F 




20225 


252.0 


62645 


5.1 


Davenport & Heintz 


4949 


23.8 


4949 


21.0 


4949 


3.8 


Hong-90 


27 


0.1 


27 


0.2 


27 


3.5 


Solotareff-3 




T/O 




T/O 


243325 


12.3 


Collision problem 




T/O 




T/O 


45979 


5.0 


Random trivariate 




T/O 




T/O 


877 


14.9 


Off-center ellipse 


4569 


153.9 


2705 


51.5 


4593 


5.0 


Concentric circles 


41 


0.2 


41 


0.2 


41 


3.7 


Non-concentric circles 


41 


0.2 


41 


0.3 


41 


3.7 


Simplified ESP 




T/O 




T/O 


56105 


10.2 


Simplified Putnum 


55021 


1026.0 


10517 


193.0 


10517 


3.9 


Simplified YangXia 




T/O 




T/O 


6313 


5.2 


Simplified SEIT 




T/O 




T/O 


F 




Cyclic-3 


381 


3.7 


381 


4.2 


381 


3.8 


Example from 


73 


0.2 


67 


0.1 


F 





Table 1: Table detailing timings and cell counts for computation of full CADs 
described in Sectional 



using the three implementations 



